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This paper considers the control of coupled aeroelastic aircraft model which is configured 
with Variable Camber Continuous Trailing Edge Flap (VCCTEF) system. The relative de- 
flection between two adjacent flaps is constrained and this actuation constraint is accounted 
for when designing an effective control law for suppressing the wing vibration. A simple 
tuned-mass damper mechanism with two attached masses is used as an example to demon- 
strate the effectiveness of vibration suppression with confined motion of tuned masses. In 
this paper, a dynamic inversion based pseudo-control hedging (PCH) and bounded control 
approach is investigated, and for illustration, it is applied to the NASA Generic Transport 
Model (GTM) configured with VCCTEF system. 

I. Introduction 

An earlier research study conducted by NASA has proven that active control of wing aeroelasticity in- 
flight can be effective in reducing aerodynamic induced drag and enhancing lift performance. The analyses 
showed that through active control of wing twist and deflection at local wing sections the overall aerodynamic 
efficiency can be improved. As a result, the novel concept of Variable Camber Continuous Trailing Edge 
Flap system is proposed . * 1 * Two sets of control actuators are employed in order to actuate the VCCTEF. 
The light-weight shaped memory ahoy (SMA) is adapted for controlling the shape of the first two chordwise 
sections of the three-section VCCTE flap. The thrid section is controlled by the electric drive motor (EDM) 
and provides the needed active wing shaping control in-flight. A preliminary static analysis has shown the 
potential efficacy of the VCCTEF system. 

This paper presents the mechanical modeling of the VCCTEF system and incorporates it into aeroelastic 
aircraft dynamics. A single degree-of-freedom vibration with tuned-mass damper (TMD) of two attached 
masses is considered. The main mass-dashpot-spring setup simulates the one degree-of-freedom vibration of 
the flexible wing, for example, first bending mode. The attached TMDs are confined in their relative motion 
and they simulate the two adjacent VCCTEFs. The TMD is well understood in the mechanical and civil 
engineering community as a device used for vibration absorber . 2,3 However, the TMD with confined motion 
is a novel concept. 

The actuation of VCCTEFs is both bounded and constrained by deflection angle and deflection rate. 
This is due to the placement of elastomer material in between the flaps to provide piecewise continuous 
motion. However, the presence of elastomer further complicates the control system design, since it renders 
confined motion between the two flaps. Previous works on bounded input and input rate, for instance 6-8 
and the references therein, offer a good framework for control designers to embark on the implementation 
aspect of the proposed approaches. A simple ’’software limiter” idea was proposed in Hess and Snell . 7 There, 
by introducing the derivative and limited integration blocks a commanded input signal can be re-generated. 
In Lin 8 a bounded input and input rate problem description was given that, in a sense, prescribes how the 
actuator should be moving. In this paper, we will utilize the bounded control approach proposed in . 6 A 
higher order state-space system was obtained by augmenting the states with the control, hence, in this setup 
the new control to the augmented system is in fact the input rate. Inevitably, the control structure in this 
case becomes a dynamic controller. 

The pseudo-control hedging or PCH 5,9 is another viable alternative control approach which is well suited 
for dealing with such practical control issue as actuator saturation. Since a dynamic inversion of the plant 
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model is involved, its application is mainly for the stable systems. However, in the context of VCCTEF 
actuation, we propose to utilize PCH modeling technique to VCCTEF model to achieve better wing shaping 
control. In this setup, the PCH based VCCTEF control can be regarded as an inner loop design that 
incorporates any mismatch from the desired control deflection and deflection rate, whereas the bounded 
input and input rate control can be considered as an outer loop design. This proposed approach is applied 
to the NASA Generic Transport Model (GTM) configured with VCCTEF system for illustration. 

II. Aeroelastic Aircraft Dynamics with VCCTEF 

Figure 1 shows the GTM with configured VCCTEF system. As shown, there are 11 flaps per wing and 
they are equally spaced, and each flap has three sections, as shown in Figure 2. Sections 1 and 2 are driven 
by shape memory alloy (SMA) and produce high lift, whereas section 3 is driven by high bandwidth electric 
drive motor (EDM) necessary to provide active wing shaping control. At nominal cruise flight, we consider 
only the section 3 is being activated. 



Figure 1. GTM with equally spaced VCCTEF 


The equations of motion for coupled aircraft rigid-body dynamics with flexible aeroelastic wing modes 
at nominal cruise flight can be described by 


X a 

X e 


Xp 



Sa 



( 1 ) 


where x a denotes the aircraft’s rigid-body states, x e consists of displacement and velocity of aeroelastic wing 
at generalized coordinates, 4 S a denotes the control surface deflection, such as elevator, and Sf denotes the 
VCCTEF 3rd segment deflection. The matrices A aa and A ee contain aircraft rigid-body and aeroelastic 
characteristics, whereas A ae and A ea correspond to aeroelastic coupling and aircraft rigid-body coupling, 
respectively. Similarly, B ae and B ea represent coupling effect between the control surface and VCCTEF. Note 
that the dimension of overall system depends on the number of rigid-body states and number of aeroelastic 
modes included in the problem setup. 
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Figure 2. VCCTEF actuation using SMA and EDM 


In order to achieve active wing shaping control for drag reduction in flight, it should be noted that 
each flap can not be operated independently. Therefore, to ensure and approximate smooth and continuous 
variation of wing shape, it is necessary to confine the relative motion between the two neighboring flaps. 
To achieve this, an elastomer material is placed between two adjacent flaps, which produces continuous 
transition for relative motion. However, presence of elastomer also introduces additional constraints on the 
performance of VCCTEF, and these constraints need to be addressed and incorporated in the control systems 
design process. The desired relative motion is often prescribed apriori and can be used for control analysis 
and synthesis. 


III. Integrating Actuator Dynamics with VCCTEF 

Figure 3 shows the schematics of VCCTEF with the elastomer modeled as a pair of spring-dashpot 
mechanism. As shown, the relative deflection between, for example, 0\ and 62 is constrained, as well as the 
relative rate between 0 \ and 62 • It is important to note that the elastomer material properties depend on 
the operation ambient conditions. In this study, we shall use notional material properties and consider them 
as potential design parameters for further analysis. 

We assume that the EDM actuator is modeled as a second order dynamics of the form 

Oi + 2 (oj n 0 i + u^Oi = LJ^Si ; i = 1, 2, • • • , m, (2) 

where Oi and Oi are output of ith actuator and represent respectively the deflection angle and rate of ith 
flap, and Si denotes the control command to the ith flap. The parameters ( and uj n are known damping 
ratio and natural frequency provided by hardware specifications. When all m flaps are interconnected with 
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Figure 3. Schematics of VCCTEF 


EDM actuators, as shown in Figure 4, the overall actuator dynamics can be described by 


k 


a + C 

-c 

0 


0i 


b + K 

-K 0 


ei 


Si 

h 


-C 

CL “h 2 C 

-C 


@2 


-K 

b + 2K -K 


e 2 


S 2 

h 

+ 

0 

-c 

a T 2C7 • • • 

y 

@3 

+ 

0 

-K b + 2K • • • 

y 

0 3 

= b 

S 3 

h 



Cf 



h 



Kf 




Scf 


(3) 

where a = 2 (uj n , b = u;^, and S c f denotes the control command to VCCTEFs. In attaining the above we 
assumed that all EDM actuator dynamics and elastomer material properties are the same. As noted earlier 
the relative deflection output from two adjacent flaps is constrained by c^, i.e. \0i-i—Oi\ < a* and \0 m \ < a m . 
In practical application, the motion of VCCTEF is confined not only by relative angular deflection, and 
because of bandwidth limitation of actuator dynamics and the presence of elastomer material, the relative 
angular deflection rate between two adjacent flaps is also bounded, i.e. | Qi-\ — 0i\ < pi and \6 m \ < p m , and 
both oti and pi are known constants. 





Figure 4. EDM driven VCCTEF 

We also assume that the actuator dynamics for control surfaces S a takes a second order form similar to 

( 2 )- 

<5a + C a S a + K a 6 a = bS ca , (4) 

where S a is the control surface deflection angle and 5 ca the control command. As noted earlier, both S a and 
5 a are bounded by motor dynamics. 

Now, combining the aeroelastic aircraft dynamics (1) with actuator dynamics (3) and (4) yields the 
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following open-loop system description, 


Xp — ApXp H- B pa 5 a H- B pe 5f 

< h + CfSf + Kf8f = bS cf / 5 x 

+ Ca^a + Ka^a = bS ca 

< Vp — CpXp 

where y p denotes the control output variables. Figure 5 illustrates the schematic of open- loop block diagram, 
which shows the aeroelastic aircraft model interconnecting with constrained VCCTEF dynamics. 


Aeroelastic Aircraft Model 


with VCCTEF 



Figure 5. Interconnection of aeroelastic aircraft model with VCCTEF actuator dynamics 


IV. A Tuned-Mass Damper Analogy 

In order to demenstrate the efficiency of VCCTEF in suppressing the virbrational motion of the wing, 
a simplified tuned mass damper (TMD) problem is considered; see Figure 6. As shown, the main mass- 
dashport-spring suite (M, ci,fci) models one vibrational mode, for instance wing’s first bending mode, and 
it is subjected to external excitation force /(£). The attached series of two small masses m are considered 
as a pair of VCCTEFs. The two dashport-spring sets (02,^2) and (03,^3) that connect the three masses 
represent the EDM dynamics and the effect of elastomers. It should be noted that, unlike the conventional 
tuned mass damper study where displacement of mass m is not constrained, here the relative displacement 
between x 2 and x% is confined, that is \x2 — xs\ is bounded. The application of tuned mass damper as a 
tuned vibration absorber is well studied and understood, however, in this study our purpose is to show the 
effectiveness of the tuned mass damper by confining relative motion between x 2 and X3. 


m 

■» | > 



Figure 6. Single degree-of-freedom constrained tuned mass damper 
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We can derive the equation of motion for the tuned mass damper problem shown in Figure 6, and it is 
given by 
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(6) 


where e = x% — X2 . 

For purpose of illustration and simulation, we assume that M = 100, c\ = 2, &q = 100, m = 2, 
C2 = C3 = 0.4, k 2 = 10, ks = 1.8, and f(t) = 10sin(£). It is clear that the natural frequency of mass M is 
1 rad/s, which also coincides with external excitation frequency. Here, we consider lightly damped system, 
though both C2 and C3, as well as k 2 and ks , can be adjusted actively by introducing appropriate control 
feedback. 

Figure 7 shows the displacement of mass M with and without the tuned mass damper, and the suppression 
of vibrational motion for mass M is apparent when tuned mass damper is activated. However, in this case 
no relative displacement constraints are imposed. Now, when relative motion between X 2 and xs is confined 
by 1, i.e., \e\ < 1, though the effectiveness of absorbing vibration motion of mass M is reduced, it is still 
acceptable, see Figure 8. And the results can be further improved by means of active control system design. 3 



Figure 7. Displacement of mass M with and without TMD 


6 of 15 


American Institute of Aeronautics and Astronautics 



Figure 8. Displacement of mass M with bounded TMD, |e| < 1. 


V. Control of Flexible Wings Using VCCTEF 

In this section, we only focus on controlling the flexible wings using VCCTEF by utilizing its model 
developed in Section III. To this end, recall the aeroelastic wing model from (1) and the VCCTEF dynamic 
model given in (3) 

f %e — A e e%e T B ee Sf /«\ 

l S f + CfSf + K f 5 f = b5 cf 

To derive (7) we have neglected aeroelastic coupling effects. It should be noted that the dimension of x e 
depends on the number of aeroelastic modes retained in the model. Note that the first equation of (7) with 
bounded Sf and Sf is precisely the problem setup considered in 6 and. 7 In this study, the actuation of Sf, and 
hence Sf, is generated by the motion of VCCTEFs, and the commanded control input S c f is so designed to 
produce desired bounded Sf and Sf. The complexity of VCCTEF dynamics suggests a better management 
and coordination of its actuation to achieve effective wing shaping control. Therefore, the control law 
architecture proposed in this section consists of two integrated loops. The deflection and deflection rate 
saturation control law developed in 6 is used as an outer loop, which is then integrated with an inner PCH- 
based controller 5,9 for managing the VCCTEF dynamics subjected to bounds and constraints. Figure 9 
shows the schematics of the proposed integrated control architecture. As shown, the control command S c f is 
a combination of bounded position/rate control and command shaping control. The command shaping for 
S c f is achieved by utilizing PCH approach 5,9 integrated with dynamic inversion, as shown in Figure 9. 
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V.A. Coordinate transformation 


Given the position and rate constraints of VCCTEF, we can define a coordinate transformation to normalize 
VCCTEF dynamics (3). First, let Sf = ESf, where E is given as 
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then, the VCCTEF dynamics can be rewritten as follows, 

Sf + C5f + ESf = BS c f , 

where C = ECfE -1 , K = EKfE -1 , and B = bE. Note that Sf is a vector of relative deflection between two 
adjacent flaps, except the ra-th flap. Alternatively, we can choose the first flap to be unconstrained. Next, 
we define the controlled outputs, Sf and Sf, as being the normalized deflection and deflection rate bounds 
by D\ and D 2 , where D\ = diag(aq, 0 L 2 , • • • , a m ) and D 2 = diag(pi, P 2 , ■ • ■ , pm )• That is, 

Sf = DpSf and 1/ = DpS f . 


Therefore, each element in Sf and Sf is being normalized by ±1. We now can rewrite Eq. (7) as follows, 

x e — A ee x e T B ee Sf 
< h + cSf + kSf = BS cf 
Sf = E~ 1 D\Sf 
K Sf = E ^Z^2^/ 


V.B. Bounded position/rate control design 

In this section, we’ll use the extended High Performance Bounded (eHPB) control algorithm proposed by 6 
to design a bounded position/rate controller for suppressing the vibrational motion of the flexible wing. The 
construction of control law is based on the standard LQR approach for augmented control system. Consider 
the aeroelastic wing model; the first equation of (8). The augmented aeroelastic wing model is derived by 
augmenting the state x e with the control input Sf, i.e. 


X = AX + BSf , 


(9) 


where 


V = 

x e 

,A = 

A-ee 

B ee 

,B = 

0 
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and I m denotes an m x m identity matrix. As shown in Section V.A, both Sf and Sf can be normalized by 
a ’’state” and ’’input” transformation. Let X = TX, where T is defined by 


In 0 

0 D^E 


and Sf = E 1 D 2 Sf. Then, (9) can be equivalently described by 

X = AX + BS f (10) 

where A = TAT -1 , B = TBE~ 1 D 2 , and each element in Sf and Sf is bounded by ±1. The problem 
described in (10) is precisely the eHPB control problem considered in, 6 and the following theorem provides 
the stabilizing control law. 
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Theorem 1 6 Consider (10). Let Q l = Q = diag(Qi,rI m ) > 0 be chosen such that (A,Q) is observable ; 
where r > 0 is a scalar. Then, the control law 


5 f = X^P(r)X ( 11 ) 

asymptotically stabilizes (10) subject to the normalized constraints, where, for a given r > 0, P(r) is the 
unique stabilizing solution to the following algebraic Riccati equation 

•4*P(r) + P(r)A ~ -P(r)BB t P(r) +Q = 0 , 
r 

and with r(X ) = min{r\X £ £(r),Wt > 0}; where 

£{r) = {X\X t P{r)X<p), 

/? = [max(diag(B t PB /r 2 ))] -1 

In summary, the control law given in (11) guarantees that the augmented closed-loop ” state” trajectories 
started within the Lyapunov ellipsoid £{r) will asymptotically converge to the origin. Though the input to 
the aeroelastic wing model is Sf, the controller derived in (11) is a dynamic one, since X consists of both x e 
and Sf. Hence, we can describe Sf more explicitly as 

S f = -K 1 (r)x e -K 2 (r)S f (12) 

where K\(r) and K 2 {r) can be extracted from (11). It is important to note that the control actuation; be it 
Sf or Sf , can not be generated instantaneously. In the case of VCCTEF, the control actuation needs to go 
through a negotiated process in order to generate a smooth flap profile for flutter suppression. 

V.C. PCH-based VCCTEF actuation 

As indicated earlier, the constrained VCCTEF dynamics in application to the flexible wing control poses a 
great challenge. Given the bounded control law (11), this section is to design an effective VCCTEF control 
management scheme to achieve desired flap profile. In this regard, the PCH-based control approach 5,9 is 
proposed, which utilizes the dynamic inversion scheme. Relatively speaking, the actuator dynamics including 
the elastomer can be modeled more accurately and its dynamics inherently stable, so we should not expect 
large mismatch from dynamics inversion process. Therefore, we consider PCH-base approach a suitable 
scheme, and it is used as an integrated part of command shaping for outer loop input command model. 

To derive the dynamic inversion for VCCTEF model, recall 


Sf + CSf KSf = B8 c f . 


(13) 


Note that B is an invertible matrix. Since both Sf and Sf are measurable, we can define 

S cf = -B-fCSf + KSf-v), 

where v is a new control input. Substituting the above into (13) yields 

5 f =v, (14) 

and we can choose v to be a simple linear position and rate controller. Because the VCCTEF motions are 
bounded and constrained, the commanded flap deflection angle S c f will differ from the actual actuator flap 
deflection angle Sf. In this case, the PCH implementation will help modify the commanded model so that Sf 
follows closely with S c f. The dynamic inversion based PCH scheme for VCCTEF model is shown in Figure 
9. 

The overview of the proposed control scheme is as follows. As noted earlier, if we discount the VCCTEF 
actuator dynamics, the bounded controller given in (11) guarantees the closed- loop stability, provided that 
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the initial conditions originate from a bounded region. As shown in Figure 9, the control law (11) serves as 
a ’’reference” generator and outputs the reference commands. The PCH signal v is defined as 


v = v - v a , 

and it is fed back as a compensation to shape the reference model. That is, 


S f = -K x {r)x e - K 2 (r)S f - K^v , 


where = & 3 D 2 and k% > 0 is a design parameter. Therefore, the commanded flap deflection and deflection 
rate are given by 

(^/)cmd — K\Ef 

(<t)cmd = A4/ 

where K\ = k\KD\ and K 2 = k 2 CD 2 , and k\ > 0, k 2 > 0 are the design parameters. 


PCH Model 



Figure 9. Dynamic inversion based PCH scheme for the control of VCCTEFs 


VI. Simulation Results 

Here, we present an application of proposed PCH-based bounded control law to the GTM aeroelastic 
wing as an illustrative example. The GTM is configured with 15 VCCTE flaps in each wing; hence, m = 15, 
and it is considered at cruise condition at Mach 0.92 and altitude of 35,000 ft. At this condition, the 
GTM has already exhibited fluttering behavior. The relative deflection angle and deflection rate between 
the two adjacent flaps are, respectively, 2 deg and 60 deg/sec. That is, |#i_i — 6i\ < 2°, \Q\$\ < 2°, and 
\6i-i — 0i | < 60 deg/sec and \0i$\ < 60 deg/sec. In addition, we consider the flap actuator dynamics at two 
natural frequencies; 10 Hz ad 30 Hz, and damping ratio of 0.8. The elastomer has the natural frequency of 
10 Hz and damping ratio of 0.8. In this study, we have chosen the design parameters fci, k 2 , and k% to be: 
ki = k 2 = 0.5 and k% = 0.1. 

For better illustration, we only consider five aeroelastic modes in this study; including the two unstable 
(anti-symmetric) wing bending modes. The objective of bounded control law design then is to suppress the 
vibrational motion of the wings using VCCTEFs. The consideration of meeting, in addition, the overall 
aircraft handling qualities/performance specifications is subject of future research. Table 1 contains the list 
of aeroelastic modal frequency and damping considered in the simulation. 

The proposed control law architecture is able to asymptotically stabilize the unstable modes when the 
VCCTEF actuator frequency is at either 10 Hz or 30 Hz. In the case of 10 Hz, Figures 10 and 11 show that 
the relative flap actuation are all bounded by ±1. Moreover, Figures 12 and 13 show the actual individual 
flap motion. Figure 14 shows the closed-loop response in the modal coordinates. In generating the above 
results we have to pay special attention to the choice of initial conditions, for the limited/bounded control 

10 of 15 


American Institute of Aeronautics and Astronautics 









Table 1. Aeroelastic modal data for GTM 


Mode ID 

Frequency (rad/sec) 

Damping 

unstable 

18.1 

-0.0025 

unstable 

37.6 

-0.0283 

bending 

20.3 

0.135 

bending 

42.2 

0.194 

torsion 

69.1 

0.0092 


authority will not render global closed-loop stability. In this regard, Figure 15 indicates the trade between 
the control effort, in terms of r, and achievable ellipsoid, in terms of ||x e ||. As can be seen that increasing the 
control effort; by reducing r, reduces the size of Lyapunov ellipsoid. In other words, large control effort can 
be implemented when closer to the ” origin,” and small gain is implemented when far away from the origin. 

One important note worth mentioning is that when we follow the conventional approach, that is, the 
bounded control law (11) is applied directly to drive the VCCTEF dynamics which in turn affects the flexible 
wing, the control law fails to stabilize the closed-loop system when the VCCTEF actuator bandwidth is lower, 
say at 10 Hz. Therefore, the advantage of PCH based control modification is the adaptation capability of 
the control system for varying actuator dynamics. 

VII. Conclusion 

In this paper, we presented a PCH-based control system design approach that utilized the VCCTEF 
model inversion integrated with the bounded input and input rate control law. Because of the complexity of 
the VCCTEF dynamics, the PCH was applied to act as a control modifier to shape the command reference 
that is achievable, given the constraints on relative flap deflection and deflection rate. A simple TMD model 
was used to demonstrate that the concept of VCCTEF is efficient in suppressing the vibrational motion, 
and the proposed control architecture can be considered, in an extent, as an active tuned-mass-damper 
for flexible aircraft systems. It was also shown that when the bandwidth of flap actuator dynamics was 
low, the conventional control approach failed to stabilize the closed-loop system, while the proposed control 
architecture was able to. On the other hand, an alternative control approach might be to take the actuator 
dynamics into account as part of plant model, and design the bounded state/input control for the augmented 
system. This is a subject of ongoing research. 
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Figure 10. Normalized relative flap deflection angle 5f time history. 



Figure 11. Normalized relative flap deflection rate 5f time history. 
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Figure 12. 


Individual flap deflection angle Sf (deg) time history. 



Figure 13. Individual flap deflection rate 5f (deg/sec) time history. 
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Figure 14. Closed-loop response in modal coordinates subject to initial conditions. 
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Figure 15. Achievable ellipsoid as function of control effort (r). 
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